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Xy^ • Monte Carlo techniques are used to investigate the equilibrium threshold con- 

centration, Xe, in the dilute anisotropic antiferromagnet FexZn\-xF2 in an applied 
Ch ' magnetic field, considered to be an ideal random-field Ising model system. Above Xe 

(y[3 ' equilibrium behavior is observed whereas below x^ metastability and domain forma- 

tion dominate. Monte Carlo results agree very well with experimental data obtained 
using this system. 
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Q ' The dilute antiferromagnet (AF) Fej.Zni_xF2 in an applied field is a realization ||l|,g] of the 

random- field Ising model (RFIM). Many studies ||^ have been done on this system for x < 0.75. In 

>" 
On ' such cases, metastability and domain formation mask the equilibrium critical behavior, particularly 

o 

ON 

o . . 

f^ ' tions. The specific heat, C„i, on the other hand, is not as greatly affected by domain formation 

o' 

-I— > 



in scattering measurements that are by nature dominated by long-range antiferromagnetic correla- 



except very close to the transition, Tc{H), since it is primarily sensitive to short range correlations 
H. Only recently has it been discovered pi that equilibrium scattering behavior can be observed 
f-H ' for X — 0.93 with no evidence of domain formation. Domain walls form with little energy cost when 

o' 

J-» ' vacancies are so numerous that magnetic bonds can be largely avoided. At high magnetic concentra- 

J^ ' tions domain walls must cut a large number of magnetic bonds and long-range order (LRO) is stable, 

C^ ' as in the Imry-Ma domain wall energy arguments [|6|. Hence, the RFIM can be studied in equilibrium 

for X = 0.93 and a transition to LRO is observed, consistent with theory ||7|. The question remains 
as to the nature of the disappearance of the domain walls as the magnetic concentration increases. 
Do they disappear gradually or is there a critical concentration above which they do not form? 
We have performed Monte Carlo (MC) simulations of the RFIM modeled as closely as possible 



after FexZni^xF2- We provide evidence that there is a critical equihbrium threshold concentration, 
Xe, above which domain formation does not occur. We estimate this concentration to be 0.75 < Xe < 
0.80. The simulation results also shed light on the experimentally observed B formation of domain 
walls upon approach to Tc{H) for x < Xe after cooling in zero field and subsequently heating in a 
field (ZFC). Finally we demonstrate that cooling in a field (PC) for x < Xe results in a metastable 
state at low T . 

The staggered magnetization [Mg] and correlation functions for the first three nearest-neighbor 
(NN) pairs have been calculated using MC simulations. The magnetic lattice corresponding to the 
body-centered-tetragonal Fei-xZnxF2 lattice is described as two cubic sub-lattices of size Lx Lx L 
each, delineated as one dimensional arrays bit coded to accommodate large lattice sizes. Most 
of the results reported here were obtained with L = 64, corresponding to more than 5.2 x lO'^ 
sites magnetically occupied with probability x. Periodic boundary conditions are used. At each 
temperature magnetic sites are visited randomly with each site visited an average of once per MC 
step. Upon each visit, a spin is flipped with a probability given by the metropolis algorithm. The 
temperature was changed in steps of 0.01 K after N MC steps. N was increased until the results of 
the simulation were largely independent of N for H — 0. Unless otherwise stated, we used N = 500 
for the results reported here. Using N ~ 1000 gave essentially identical results for H = 0. The 
temperature scanning procedures correspond to the experimental ones, i.e. cooling and heating for 
H = 0, and PC and ZPC for H = IS T. Fei-xZnxF2 is well represented |^ by the Hamiltonian 

H=J2j,,S,-S,~Dj2sf-H-J2s. . (1) 

In the simulations we use the Hamiltonian 

H=J2 J^JS^'S/-HY,S^' , (2) 

which corresponds to the Ising limit D —> oo. The first three NN exchange interaction strengths 
are taken from spin-wave dispersion measurements [p[ on FeF2: Ji — —0.069 K; J2 — 5.278 K; and 
J3 — 0.279 K. All other interactions are negligible and are not included in the simulations. The 
transition temperature increases with anisotropy in this system pO[. Hence, the infinite anisotropy 



of the model results in a transition temperature much higher than in the real system. As is done 
elsewhere [ni|, the included interactions were all scaled by a factor 2/3 to make the transition 
temperatures more closely correspond to the ones observed [n2| in Fei^xZnxF2. 

Figure 1 and its inset show typical behavior for x — 0.60 < x^ for three cases: H = 0, which 
is free of hysteresis; FC with H = 13 T; and ZFC with H = 13 T. In the inset, the behavior 
of Cm = dE/dT vs. T is shown. For H = random-exchange behavior is observed with slight 
rounding from finite sample size effects. Taking the rounding into consideration [[l3| , we estimate the 
transition temperature Tn ~ 45.0 K. For the ZFC and FC procedures, the specific heat is depressed 
in T as expected. Away from the shifted transition, the two behaviors are essentially identical with 
a shape more symmetric than for ff = 0, as is observed in experiments B. Since the ZFC and 
FC behaviors are the same outside the distorted regions, they clearly share the same Tc{H). On 
a plot of the specific heat versus the logarithm of the reduced temperature, Tc{H) appears close 
to the peak of the ZFC data, w 40.5 K. From the main figure, we see that for iJ = 0, Ms falls to 
zero very near the Tn taken from Cm- Upon FC, however, Ms remains unusually small quite far 
below the Tc{H) determined from Cm- Upon further cooling, the system begins to rapidly evolve 
toward AF LRO with a curvature opposite to that for iJ = 0. Finally, at a rather well defined 
temperature (Ti) the Ms curve abruptly changes curvature. Upon further cooling, Ms increases 
smoothly and slowly, much like the H — case but shifted to lower T. Upon ZFC, Ms vs. T follows 
a much smoother evolution but nevertheless approaches Tc{H) with a curvature opposite to the 
H = case. The unusually rapid disappearance of AF LRO well below Tc{H) is attributed to the 
formation of domains and has been observed experimentally S . The FC behavior in experimental 
systems, on the other hand, is quite different from that shown in Fig. 1. In the experiments, the 
domain formation prevents significant Ms down to low T . Hence, we examined this behavior further 
using the MC techniques. We found that by slowing the temperature scans by using N = 1000 with 
L — 64, the domains that form for x < Xe remain stable to much lower temperatures. We also did 
a simulation with i = 128 and N — 500 and again the domains were stable to lower temperatures. 
Extrapolating these two trends to the much larger experimental systems that are cooled much more 



slowly relative to the spin flipping time, it is not surprising that the experimental systems show 
metastable domains that are very stable at low temperatures. 

Figure 2 shows the ZFC and FC behavior for x = 0.60, 0.70, 0.80 and 0.90 at i/ = 13 T. It is clear 
that the hysteresis decreases as x increases and is absent for the cases x > 0.80. It thus appears 
that the systems at larger x are in equilibrium. To demonstrate that there is a critical threshold 
concentration above which equilibrium behavior, such as that observed for x > 0.80, prevails, we 
plot in Fig. 3 the temperature AT — T2{H,x) — Ti{H,x), where T2{H,x) is the temperature at 
which the ZFC joins the FC one just below Tc{H). Note that the ZFC and FC curves separate at 
Ti{H) < T2{H). It is clear that AT becomes zero somewhat above x = 0.70 and below x — 0.80, 
strongly indicating a critical concentration 0.75 < Xe < 0.80, which is much larger than the magnetic 
percolation threshold Xp — 0.25. The field dependence of Xe is weak. Within the accuracy of the 
present MC resuhs, we find that for iJ = 7 T and i? == 13 T, 0.75 < Xe < 0.8. This suggests that 
the critical concentration is primarily geometric in origin. Domain formation in the ZFC procedure 
indicates that for iJ > 0, the free energy must favor domain formation close to Tc{H) for x < Xe- 
Such behavior is indicated in local-mean-field simulations J14| . However, there is no evidence for a 
latent heat associated with the temperature at which the free energy difference for the domain state 
and AF state switches sign. 

Finally, we demonstrate in the inset of Fig. 3 that the ZFC peak in d{Ms'^)/dT (the Bragg 
scattering) vs. T for x ~ 0.60 is at a lower temperature by more than 1 K than Tc{H) as determined 
from the Cm- This behavior is in agreement with neutron scattering results H but in contradiction 
to the so-called 'trompe I'oeil' phenomenology Il3]. 

In conclusion, we have presented MC evidence for the existence of a critical concentration x^ above 
which equilibrium RFIM behavior prevails and below which domains form close to the transition. 
The numerous vacancies for x < Xe allow domain formation without a great energy cost. We note 
that the vacancies themselves percolate at a; = 0.75 and surely Xe is related to this percolation 
threshold, although this may not be the entire story. Domain walls take advantage of surfaces 
of vacancies, not just the filiamentary percolation threshold structures. Hence, it would not be 



surprising if Xe were somewhat different from the vacancy percolation threshold. Further simulations 
are currently under way at i = 128 and should help to determine Xe more precisely as well as 
determine the influence of the finite sample size on the present results. Certainly, critical behavior 
studies should be performed in the real systems only for x > Xe- 
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FIG. 1. Ms vs. T for a: = 0.60 a.t H = 0, and for ZFC and PC at H = 13 T. The transition temperatures 
are determined from the specific heat behaviors shown in the inset where the units are expressed in terms 
of two magnetic sites per unit cell. In both the main figure and inset, curves are, from left to right, for the 
FC, ZFC and _ff = procedures. 
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FIG. 2. The ZFC and PC Ms vs. T ior x = 0.60 (a), 0.70 (b), 0.80 (c) and 0.90 (d). The upper curves 
are ZFC in each case except x = 0.80 and 0.90, where the ZFC and FC curves coincide at all temperatures. 
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FIG. 3. The temperature range over which hysteresis occurs ai H — 13 T, AT = T2{H,x) — Ti{H,x) 
vs. X. The inset shows d{Ms'^)/dT (arbitrary units) and Cm/R vs. T for x — 0.60. Note that the peak in 
d{Ms'^)/dT vs. T is significantly lower than the corresponding peak in Cm/R- 



